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In Part I (physics/0405091), it was shown that it is possible to achieve efficient heating of cylindri- 
cal targets by electromagnetic beams with transverse dimensions that are several orders of magnitude 
greater than those of the cylinder. 

Part II contains derivation of the detailed conditions of efficient heating in the longitudinal ge- 
ometry and establishes a broader domain of parameters providing efficient heating in the transverse 
geometry. 

One possible implementation using currently available technology is considered: a nanotube 
(length 1 cm, diameter 15 nm) may be heated in the transverse geometry by a 30 femtosecond 
pulse of a 50 GW Ti:Sapphire laser (wavelength 800 nm) to a temperature of 1 keV with heating 
efficiency of 16 per cent. The resulting energy density exceeds 100 MJ/cm3. 

An exact formula in the form of a one-dimensional integral is obtained for the amount of energy 
absorbed in the cylinder heated by a gaussian electromagnetic beam in the longitudinal geometry. 
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I. INTRODUCTION 

In Part I (Ref. [l]), it was shown that it is possible to 
achieve efficient heating of cylindrical targets by electro- 
magnetic beams with transverse dimensions that are sev- 
eral orders of magnitude greater than those of the cylin- 
der. 

Part II contains derivation of the detailed conditions of 
efficient heating in the longitudinal geometry and estab- 
lishes a broader domain of parameters providing efficient 
heating in the transverse geometry. 

One possible implementation using currently available 
technology is considered: a nanotube (length 1 cm, di- 
ameter 15 nm) may be heated in the transverse geometry 
by a 30 femtosecond pulse of a 50 GW Ti:Sapphire laser 
(wavelength 800 nm) to a temperature of 1 keV with 
heating efficiency of 16 per cent. The resulting energy 
density exceeds 100 MJ/cm3. 

An exact formula in the form of a one-dimensional in- 
tegral is obtained for the amount of energy absorbed in 
the cylinder heated by a gaussian electromagnetic beam 
in the longitudinal geometry. 



II. DOMAINS OF PARAMETERS PROVIDING 
HIGH ABSORPTION EFFICIENCY IN THE 
LONGITUDINAL GEOMETRY 



In Part I of this work (Ref. the following asymp- 
totic formula for absorption efficiency in the longitudinal 



geometry was derived: 
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Here rj is the absorption efficiency, W a is the power ab- 
sorbed in the cylinder, W is the power in the incident 
beam, a is the radius of the cylinder (a system of units 
is used where the wave vector in free space ko = — = 1), 
Pi ~ 77 "C 1, T\ is the beam waist radius defined by 
e-fold field intensity attenuation, 
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e is the complex permittivity of the cylinder. One can 
see that the refracted fields and, therefore, the absorption 
efficiency do not change when p2 is substituted with — pi- 
Therefore we shall assume that 
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(p2 can only be real if e is real and, therefore, there is 
no absorption in the cylinder, so we do not consider this 
case) . 

It should also be noted that within the accuracy of the 
approximations used to derive Eq. (TTJ) it is possible to 
substitute P2 in that equation with \JJs — l)a. 

In this section, we are going to find the domains of 
parameters that provide high absorption efficiency in the 
asymptotic case p\ <C 1, I — \ln{p%)\ ^> 1. We consider 
absorption efficiency "high" if it is at least of the same 
order of magnitude as j . Although the value of the latter 
expression tends to zero as p\ tends to zero, it decreases 
very slowly as r\ increases and is quite high (typically 
tens of percent) for any realistic values of t\. We are 
only interested in parameters that provide high efficiency 
when p\ varies by an order of magnitude, as if absorption 
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efficiency is oniy high for one value of p\ , it may be low 
for the relevant Gaussian beam, which is a linear com- 
bination of cylindrical waves with different values of p\ 
(Eq. ([T]) was derived using the one- wave approximation) . 

The right-hand side of Eq. ([!} strongly depends on 
three parameters: the real and imaginary parts of com- 
plex permittivity e and the radius of the cylinder a. Anal- 
ysis of Eq. {!]) may be broken down into several cases 
depending on the value of P2 ■ Essentially, depending on 
this value, we use four different approximations for the 
Bessel functions: for p% < 1, for p' 2 ' > 1, for p2 in the 
vicinity of a positive zero of function J\{x), and, finally, 
for p 2 in the vicinity of a positive zero of function Jq(x). 
Evidently, the domains of applicability of the relevant 
approximations partially overlap. 

This analysis is rather cumbersome, but it has some 
elements common to all four major cases. These elements 
may be illustrated as follows. If 



calculated using the approximations will provide high ab- 
sorption efficiency. On the other hand, this lack of accu- 
racy does affect the essential formula for the second case: 
\X' — 1\ < 1, so high absorption efficiency will be provided 
only in the upper part of the wide range of Eq. ([6]), if the 
optimal values are calculated using the approximations. 
This defect of the analysis is somewhat mitigated by the 
fact that / cannot realistically be much greater than 20. 
One may also question the value of asymptotic formulae 
derived using the approximation I 3> 1 (that means that 
the Gaussian beam is several orders of magnitude wider 
than the cylinder). The author believes, however, that 
these formulae will be a very useful guide to more real- 
istic cases. In particular, it should be emphasized that 
for smaller I it is generally easier to achieve high absorp- 
tion efficiency, so the formulae show that the longitudinal 
geometry may have great practical value. 
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For the purposes of illustration, let us treat X' and X" 
as independent parameters. The standard procedure of 
maximization over, say, X" is complicated by the fact 
that I is defined to an accuracy of unity (in the expansion 
of the fields of a Gaussian beam into cylindrical waves, 
values of pi differing by an order of magnitude make a 
significant contribution to the power of the incident beam 
- see Eq.(22) of the first part of this work (Ref. [![)). As 
a result, we have to consider two cases: \X' — l\ 3> 1 
and \X' - l\ < 1. If \X' - l\ > 1, then X" = \X' - l\, 
in the maximum, and the maximum is rather 



broad (its width at the level of 50% of the maximum 
corresponds to variation of the argument by at least an 
order of magnitude). Obviously, rj > i if \X'\ < I. On 
the other hand, if \X' - l\ < 1, then X" - 1 and 77 - 1 
in the maximum. Moreover, rj > 7 if 



< X" < I. 
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It is important to emphasize the difference between these 
two cases from the point of view of the subsequent anal- 
ysis. For the analysis to cover all possible values of P2, 
the approximations for the Bessel functions are used in 
relatively broad ranges. Therefore, for some values of 
P2 the worst accuracy of the approximations for func- 
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jj^t may be, say, 30%. Nevertheless, the essential 
formulae for the first case (e.g. X" - \X' - l\, \X'\ < I, 
rj > j) still hold true in the sense that the optimal values 
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thus 
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As we assume that p\ <C \P2\-, Eq. © yields 
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therefore, 
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Let us define real values x, y, and e" by the following 
equalities: 



and 
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Then 
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As pi is denned with an "accuracy" of an order of mag- 
nitude, it may be roughly said that I is defined with an 
accuracy of unity. Let us determine to what accuracy 8 



value -^r is defined: 
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so 
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Eq. (T3J) yields 



4 ;2 2 £ " 
77 - -Z 



In the maximum 
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and 
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Condition 77 > j is met if < Z 2 . On the other hand, 
|p 2 | < 1, or |(e - l)a 2 | < 1 (see Eq. ©), thus 
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Condition 77 > j is met if e" > 1. Thus, we obtain the 
following conditions of high absorption efficiency in this 
case: 
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On the other hand, 
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2. |a;-l|»i 

In this case, absorption efficiency (see Eq. ([13])) has a 
maximum when 



4/2 



a 4 Z 



(x-l) 2 (l-a 2 Z) 2 ' 



(17) 



Thus, if x w 1, then 1 a 2 Z, Z -4-, value 1 — a 2 Z is 
defined with an accuracy of a 2 AZ ~ a 2 ~ j, x is defined 
with an accuracy of -371 ~ j. Therefore, if |x — 1| 3> j, 
then value is sufficiently well defined, and value 
is defined with an accuracy of 
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(see Eqs. l|T5lT7|l ). 

If, on the other hand, \x — 1| ~ j, then 
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(we have taken into account Eq. ([T5f ). 
Condition 77 > 7 may only be met if 
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Let us consider the following cases. 

a. x < 1. If x < 1, Eq. (j2"6"|) is equivalent to the 
following equation: 
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In fact, if Eq. (|2T|) is true, then 
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(as x w 1). Let us consider this case. 
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If, on the other hand, Eq. (127)) is false, and \y\ ^> (a-i) 2 > 
then 
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as x < 1 and x > 0. 

Thus, the conditions of high absorption efficiency are: 
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plus Eq. ([ST]), and 
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(see Eq. ((9])), we have \y\ < 4j, e" < 4j, and, in view of 

Eq. m < 1. 

Taking into account Eq. (|30p . we obtain 
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Thus, if \y\ < then, in view of Eq. (l32|) . condition 
|jP2 1 < 1 is satisfied. 

b. i» 1. In this case we have Eq. (|3"2"j) again. As 
i» 1, we also obtain 
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From Eq. (|26[) we have 
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Eq. (|32|) may be rewritten in this case as 
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Eqs. (|42I43|) imply \p 2 \ -C 1. In fact, x > 1, therefore 
4j 3> /; on the other hand, from Eq. (|4"2")1 we obtain 
\y\ ~ 1, so |y| < 4 J . On the other hand, |^-| - 1, 
^ < 1, so Eq. J43|) yields e" « 4,. 

It is also necessary to check that in the expression for 
D (Eq. (136) of Part I of this work (Ref. Q)), the term 
with \ is at least of the same order of magnitude as the 

term with p\ . 



The coefficient of 4r in D equals 
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1 J'M 

P2 Jl{V2) 



h = -i + h 



k + i|<- 

X 



(44) 



(45) 



(46) 



(47) 



where x = -4™. We may express e in the following form: 



e = — 1 + a 2 / a , 
where /<* = Za', la'l < 1. Then 



Eq. © yields 
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Let us obtain an approximation for / 2 , taking into ac- 1. p" > 1 

count that |p 2 | "C 1: 

The approximation of the Bessel functions for an ar- 
1 J[(p 2 ) 1 | gument with a large magnitude may be used: 
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We have taken into account that, e.g., substitution of 
J[{p2) by \ introduces an error equivalent to introduction 
of a factor (1 + a"pi)i where a" ~ 1 and \p 2 \ <C 1. 
In view of Eq. (|46"]) . the coefficient of ^ in D equals: 
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This coefficient is typically of the same order of magni- 
tude as I (~ I). We may assume that it is > 1, as I is 
defined with an accuracy of unity. 
The coefficient of p 2 m D equals 
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Therefore, a sufficient condition ensuring that the term 

with -4j in D is at least of the same order of magnitude 

p 2 

as the term with p\ may be written as follows: 



MP2) _ cos(p 2 - f ) r sin(p 2 - if) = 

Jl(P2) ~ COs(p 2 - ^) ~ COs(p 2 - x) 

/ 3tt\ tan(z) + tan(zp") 
= - tan r _ Tj = -l-tan(z)tan( lp ")' (M) 
where z and p" are real, and 
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tan(ip") : 
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Let p' = Re(p2), so P2 = p' + ip" ■ Then 
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thus, 
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Let us consider the following case: 
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Let us prove that 
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In fact, if \p'\ > p", then 
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on the other hand, if \p'\ < p", then 
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on the other hand, if \p'\ < p" , then 
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Therefore, in Eq. (|66p . we may neglect value //s in 
comparison with |ln(pi)|, as the latter is defined with 
an accuracy of unity. We may also neglect value , 2 + p „2 
there for the following reason. Let us assume that 
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we assume that this is too little (we are interested in 
values rj > j). 

While the magnitude of the real part of the denomina- 
tor in Eq. (|6"6")) may be much smaller than I, that would 
imply 
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so in this case we also may neglect ,,i in Eq. 
Therefore, we obtain 
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Let us introduce the following notation: a — [3 
Again, I = |ln(p x )|. As p" > 1, 
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In view of Eq. @, we obtain 
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Evidently, a/3 > 0, otherwise e" < 0, and there is no ab- 
sorption. From the definitions of a and (3 and Eq. ([3]) we 
then obtain a > and f3 > 0. High absorption efficiency 
may only be achieved if 
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whereas if a < [3 ,then 
Similarly, 
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As I is defined with an accuracy of unity, let us consider 
the following two cases (as I cannot realistically be much 
greater than 20, we shall commit a sin against mathemat- 
ics and assume for the sake of simplicity that if \(3— 1\ <C I, 
then \P - l\ < 1). 

a. |/3 — Z| ~ 1 In this case 



< a < I, 



(87) 



as otherwise rj <C j . We may assume that \ < a < 1 (the 
combinations of a and /3 that we thus omit (|/3 — /| ~ 1, 
a r~i I s ) are considered together with the next case). 
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Let us introduce the following notation: 

u = a 2 - ft 2 , 



Then we obtain 
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(89) 



(90) 

as ft = I + A/3, where Ap < 1, so ft 2 = I 2 + 2lAp + A|. 
We also obtain 

2 < v < 21. (91) 
As p" > 1, then ft = £j > , so 



I > 



2a 2 ' 



(92) 



In view of the results for the next case (see below), we 
may write 2 < v < 2l 2 . Thus we obtain the following 
conditions for y and e" (cf. Eqs. (j8lM8T I59l) ) : 



\y + 4a 2 l 2 \ < 4a 2 



(93) 



4Z<e"<8a 2 Z 2 . (101) 
2. Parameter p2 in the vicinity of a zero of Ji (x) 

Let us now consider the case where P2 — Po + Ap, po 
being a nontrivial zero of function Ji(x), so po is real 
and positive (negative zeros are not relevant in view of 
Eqs. (|2"15|> ). and |Ap| < 1. For all x we have 

4(x) = -J 1 (x), (102) 



Jo(x) = ^ + J[{x), (103) 



h{x) = J^x) = ^M_±W + j»( x) (104) 



SO 
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b. \ft — l\ ~ Z. In this case 
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as otherwise rj ~ |£ <§C j (cf. Eq. ([85]) ). We may assume 
that in this case a ~ Z and /? < Z, as the combinations 
of a and ft that we thus add (\ft — l\ ~ 1, a ~ Z) were 
omitted in the previous case. Let us further divide this 
case into the following two cases: 



a ~ I, ft ~ Z. In view of Eqs. 



we obtain: 
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conditions for ?/ and e" : 
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(we left more terms in the expansion of the denominator 
as the first term there is zero). 

Let us rewrite Eq. fTJ) in the following form: 
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a ~ Z, ft <C Z. In view of Eq. ([75]) . we obtain u ~ Z 2 , 
4j- < ii -c 2Z 2 . In view of the conditions for the previous 



case, we may write for this case: u 
We obtain for y and s": 



y~4a 2 l 2 , 
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(again, Z 
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(109) 



(110) 



8 



r\ i=s Im- 



1 

A /J 



1 Ap \ I 1 , p 



2 po y VPo 2q2 
= Im — 



Im- 



1 1 i 14 



Ap ^po ^ la? ) t -r 4a 2 T 2a 2 

4 . 4Ap 



1/1 I Po 



= Im 



, 1 + 3 _ I 

Ap \p 1 2a 2 J 2p* ' iaJ 

4Ap(Z x + (Z 2 - l)Ap*) 



= Im 



= Im 



^Zi + Z 2 - J Zi + (Z 2 - Z)Ap 



4Ap"Zi 



(Zj + (Z 2 - 0Ap)(Z x + (Z 2 - /)Ap*) (Zj + (Z 2 - Z)Ap') 2 + (Z 2 - Z) 2 Ap" 2 

4Zi Ap" 



(^2-0 2 



an) 



where Zi, Z 2 , Ap', Ap" are real, 

7 1 , PO 

^1 = 



Zo 



Po 2a 2 ' 

2 2p 2 + 4a 2 ' 
Ap' + iAp" = Ap. 



(112) 



(113) 



(114) 



As Z is defined with an accuracy of unity, let us consider 
two cases: 

a. | Z 2 — l\ ~ 1. In this case we obtain 



Zi >/, 

4ZiAp" 
' Z\ + Ap"- 



(115) 
(116) 

(117) 



If Ap" < 1, then jy < -4. < i. If Ap" 
Eqs. (111211151) yield Z x « §Zp , so 77 ~ i if Ap" - 1, 
Po - 1, |Ap'| < 1. 

b. |Z 2 — l\ 3> 1. As in this case 



2p 2 4a 2 



- l\ » 1, 



we obtain 



Let us consider two cases: 

-tt > Z. In this case we obtain 



(118) 



(119) 



(120) 



4a 2 ' 



Zi 



Z 2 -Z 



> 



Po, 



Zi 



Zo-2 



Ap' 



> 



7 Po 



p a 



(121) 



(122) 



(123) 



(124) 



(we use reasoning similar to that after Eq. |86|) and 
assume that |Z 2 — l\ ^> 1 implies |Z 2 — l\ > I, so 
|Z 2 -Z|~Z 2 ). 
If Ap" < 1, then 



, 2 1 4 « ^ 4 
77<4p a 2 — = < - 

Po Po I 



(125) 



1, then If, on the other hand, Ap" ~ 1 , then 



4a 2 4 
Po ~ I ' 



So the absorption efficiency is high if 
Po I ' 



therefore, 



and 



a 2 Z~ 1 



Po ~ 1- 



(126) 



(127) 



(128) 



(129) 
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4y <c /. In this case 



and 



4Zi 



\Z 2 \ «/ 
Ap" 



Let us determine the accuracy with which expression 



(130) 
(131) 



is defined: 



Zx 



Zo-l 



Zi 



Zi 



Z 2 -IJ Z2-1-1 Z2-1 

Zx Zx 

~ (z 2 - ly ~ p ■ 

Therefore let us consider two cases: 



Zi 



In this case 



In the maximum 



4Zx Ap" 



I 2 (Zx' 



Ap" 



4Zx 1 

V ~ — — — ~ 1. 



There should also be 
\Ap'\ 

Condition 



I 2 2f 



Zx 



Zo-l 



< 1. 



Zx 



Zn-l 



< 1 



is equivalent to the following: 



0_ Po 



I 2a 2 / 



< 1. 



Ap' 



> 



; 2 • 



In the maximum we have 
Zx 



Ap" 



Zn-l 



AZx 



Ap' 



» 



Zx 



1 



\z 2 -i\ 



1 + Ap'^d 



2 1 

7 



1 + Ap' 



(132) 



(133) 



(134) 

(135) 
(136) 



(138) 



(139) 



(140) 



(141) 



In view of Eq. (|140[) we may assume that 
Ap >-, 



so there should be 



Po 
2a 2 



<l. 



(142) 



(143) 



(144) 



In view of Eq. (fHTj) we have \Ap'\-± < 1, or |Ap'|^ < 
1, otherwise -q <C j. 

3. Parameter P2 in the vicinity of a zero of Jo(x) 

Let us consider the case where P2 = go + Aq, qo being a 
zero of function Jo(x), so qo is real and positive (negative 
zeros are not relevant in view of Eqs. (|2Jt3j) ) , |Aq| < 1, 
Aq = Aq' + iAq", where Aq' and Aq" are real, and 
Aq" > 0. Using Eqs. (|102til03jl . we obtain 



Jo(qo) = 0, 

J o(<?o) = -Jx(qo)- 

Therefore, 

J (q + Aq) J Q (q ) + AqJ^(q ) 



(145) 



(146) 



Ji(q + Aq) 



■A (go) 



1 P2 ^ 1 go 
p 2 2a 2 q 2a 2 



-Aq, (147) 



(148) 



(137) We leave fewer terms in the expansions than in some of 
the above cases, as function j4^j , obviously, varies more 

slowly in the vicinity of zeros of function Jo(x) than in 
the vicinity of zeros of function Jx (x) . 



Let us denote Yi 



- + then 

qo 2a 2 ' 



rj « Im 



Jo(P2) 

Jl(p2) 



( J_ 4. P2_ 

\pi T 2a 2 



Im- 



= Im- 



4((-Aq' + iAq")Yx - I) 



Im- 



-Aq'Yx - I) 2 + (-Aq"Yx) 2 
4Aq"Fi 



(149) 



(Aq'Yx + I) 2 + Aq" 2 Y 2 ' 
There should be Y\ ^> 1, as otherwise r\ < h. Therefore, 



9o 

qo ' 2o 2 
Let us consider two cases: 



go 
2a 2 



» 1. 



(150) 
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a. \A.q' + l\ ~ 1. In the maximum 
A<? 2^ ~ X ' 



»7 ~ 1, 



(e — l)a 2 , ^2 = p' + ip" , p' and 



As |Ag'| < 1, 



go 
2a 2 



>.l. 



(151) 



(152) 



(153) 



(154) 



(155) 



b. |Ag' J% + Z| 3> 1. In the vicinity of the maximum 



p\ = (e-l)a 2 +p 2 
p" are real, p" > 0, 

e = e' + ie" is the complex permittivity of the cylinder, 
e' and e" are real, 

I = 4- 

2/ = £'-l, 

Po is a positive zero of the Bessel function J\ (x) , 
Ap = Ap' + iAp" — p 2 — po, Ap' and Ap" are real, 
Z x = - 1 1 po 



Z, = 



1 2a 2 ' 
1,3 



2p 2 1 4a 2 ' 

go is a positive zero of the Bessel function Jo (a;), 
Ag = Aq' + iAq" = p 2 — go, Aq' and Ag" are real, 
Y\ — — + 7?% • Let us now summarize the conditions of 

1 qo 2a 2 

high absorption efficiency using a hierarchical structure: 
!• H < 1 



1.1. 



1.2. |ar-l| > 
1.2.1. 



1 



1 

7' 

e" ~ /, 



(163) 



|Ag'^ + /|' 



(157) 



<1, 



|A^|<<, 



as otherwise 



T] < 



We may assume that 



(158) 



(159) 



(160) 



as Ag"^ > 1. Therefore, |Ag'| < Ag". On the other- 
hand, Eq. (|156p yields 



2a 2 l 
go 



< 1. 



(161) 



(162) 



C. Domains of parameters providing high 
absorption efficiency (summary) 

For easy reference, let us summarize the notation: 

a is the radius of the cylinder (a system of units is used 
where the wave vector in free space fco = — = 1), 

Pi 1=3 77 *C 1, where r\ is the beam waist radius defined 
by e-fold field intensity attenuation, 

I = |lnpi| < 1, 



1.2.2. 



max 



2- \P2\ > 1 
2.1. p" > 1 
2.1.1. 



2.c 



x — 1 



1 2/1 



|:r- 1| » y, 

2cc 

' (i- 1) 2 Z 

< x 



2x 



x — 1 



(x-iy 



\y\ < - 



a; > 1, 

\y + 2\<-, 

x 

2x 



{x-l) 2 l 
1 



1 > 1 



\y + 4a 2 l 2 \ <4a 2 /, 

8a 2 < e" < 8a 2 Z 2 , 
1 



I > 



2a 2 



2.1.2.1. 



Ivl < 4a 2 Z 2 , 
e" ~ 8a 2 l 2 , 

~ 2a 2 • 



(164) 



3 £ -T- ( 165 ) 



(166) 



(167) 
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2.1.2.2. 



y~4a 2 / 2 , 
4/< £ "<8a 2 / 2 . 



2.2. p 2 = Po+Ap, Ji(po) = 
2.2.1. \Z 2 -l\ ~ 1 



)) = 0,po>0, Ap = 



(168) 
Ap'+iAp" 



4a 2 ' 




-1, 




Ap' 




-1, 


Po " 


-1, 



2.2.2. \Z 2 -l\ > 1 
2.2.2.1. 4r >l 



3 

I^ 2 



2.2.2.2. 4, < 
2.2.2.2.1 



» 1, 

Ap" ~ 1, 
Ap' < 1, 
a 2 l ~ 1, 
Po ~ 1- 



^I + AP' 



4a 2 Z 



»1, 



^7 + ^ 



a^ 

Z 2 ' 

Ap'^f, 
2a 2 / ~ ' 



2.2.2.2.: 



^I + AP' 



»1, 



Ap". 



Zi 
>-p, 



2a 2 ~ ' 

. . ,,2a 2 Z 

Ap 

Po 



<1- 



2.3. p 2 = go + Ag, J (go) = 0, g > 0, Ag = 



(169) 



(170) 



(171) 



(172) 
Ag'+iAg" 



2.3.1. |Ag'^ + 



l\ ~ 1 



Ag 



1, 
1, 



2a 2 
go 
2a 2 
|Ag'| < 1, 
Ag" < 1. 



(173) 



2.3.2. |Ag'^ + l\ > 1 



Ag" 



Ag' 



»1, 
2a 2 / 



9o 



go 
2a 2 



(174) 



III. A RELATIVELY SIMPLE EXACT 
FORMULA FOR POWER ABSORBED IN THE 
CYLINDER 

Rather surprisingly, it is possible to derive an exact 
formula expressing the power absorbed in the cylinder as 
a one-dimensional integral. It should be noted that this 
formula has an obvious drawback: it gives an expression 
for - power absorbed in the whole infinite cylinder 
(— oo < z < oo). On the other hand, as we are inter- 
ested in high specific power input, it is more important 
for practical purposes to know the power absorbed over 
a finite section of the cylinder, e.g. the section where the 
linear absorbed power (absorbed power per unit length) 
is at least half as large as in the maximum (Wq 5 ). It 
seems, however, that this formula may be quite useful 
for testing of more detailed calculations. Furthermore, 
it is possible to get a rough estimate of ratio ^fe 6 -- It is 
well-known that z-dependence of the magnitude of the 
Poynting vector on the axis in a Gaussian beam is de- 
scribed by factor w a 1 +z2 ( as always, we use the system of 
units where the magnitude of the wave vector in vacuum 
equals unity). Thus, this factor is half as large as in the 
maximum at z = ±Wq. Therefore, it is reasonable to 
believe that 



J — w 



dz 



OO Wq-\-Z Z 



(175) 



Direct calculations typically yield values around 0.6 for 
this ratio. 

So let us derive the formula. Again, we describe the 
incident Gaussian beam by z-components of the electric 
and magnetic Hertz vectors: 



IL(p,<p,z) 



[ ah(X)J 1 {Xg) 
Jo 



dX, 



(176) 
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n'{p,<p,z)= [ 

Jo 



1 

l3h(X)J 1 {Xg)dX, (177) v 2 (p,(p,z)= I b 2 {X)h{X)J 1 {X 2 g) dX, (180) 



where a = 1, (3 = —i, 

X 2 

h(X) = w 2 exp(i(<p + jz)) exp(-u^y ), (178) 



and 7 = VI- A 2 (sec Eqs.(22) of Part I of this work where and ^ are giyen by Eqs .( 99; ioi) f 

(Ref. Ul)). Then z-components of the electric and mag- ( Ref fj[) for n = 1, e t = p x = p 2 = 1, e 2 = e, Pl = Xa, 

netic Hertz vectors for the refracted field are: and A 2 = e - 7 2 , p 2 = X 2 a. These Hertz vectors yield 

,-i the following refracted electric and magnetic fields on the 

u 2 {p 1 ip,z)= J a 2 (X)h(X)Ji(X 2 g) d\ 7 (179) surface of the cylinder (p — a): 

I 



H(a,<p,z) = J a 2 (X)h(X)^J 1 (p 2 ),iX 2 J[(p 2 ),(^ dX, (181) 



E{a,<p,z)= [ a 2 (X)h(X)\ l -^J' 1 (p 2 ),-^J 1 ( P2 )Aji(P2)\ dX, (182) 
Jo I £ ea e ) 

H'(a,<p,z) = J b 2 (X)h(X){i 1 X 2 J[( P2 ),-^J 1 { P2 ),X 2 2 J 1 ( P2 )} dX, (183) 

E'(a,ip,x) = J 6 2 (A)/i(A)|-iji(p2),-a 2 J((]?2),o| dX (184) 

I 



(cf. Eqs. (78-79) of (Ref. [l|). The total refracted electric The total power absorbed in the cylinder equals 
and magnetic fields are E rfr = E + E' and H r f r = 
H + H' , correspondingly. The averaged p-component of 
the Poynting vector at the surface of the cylinder equals 



f27T roo 

WS a = -a I d V dz 



n 



i^Re [^(a, V , ,) x n v , ,)], . ( W - W ) ■ P») 

= lj- Re ( E v rH ^ r *- E z fr H^ fr *). (185) We obtain: 
I 



E; fr = J {-a 2 (X)^J 1 (p 2 )-ib 2 (X)X 2 J{(p 2 ))h(X)dX, (187) 

E r J r = [ a 2 {X)^-J 1 {p 2 )h(X)dX, (188) 
Jo £ 

H r J r *= f {-ia* 2 (X)Xy[*{p 2 )-b* 2 (X)ljl(p 2 ))h*(X)dX, (189) 



fff"*=/ 6;(A)A; 2 Jr(p 2 )/i*(A)dA. (190) 
I 







We also have Integration with respect to ip in Eq. (|186p yields factor 

el 



,dX=- [°...^= (191) 



/o 



A ./n A 
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2ir. The resulting three-dimensional integral may be re- soning: 
duced to a one-dimensional one using the following rea- 



dz I 

—oc JO 



r-1 /> 1 /> 1 /> OO 

511(7) exp(i7z) c?7 / .g 2 (7') cx P(H7' z ) c h / = / / d^di gi{~i)g 2 {l') \ dzexp(i(j - ^)z) = 

JO JO Jo J-00 



10 Jo 



dldi '51(7)52(7') ' 2tt5(7 - 7') = 2tt / ^751(7)52(7). (192) 

/n 



Therefore, 



w 2 exp ( -Wo 4") 7 



X (A)d7 = -ay W ^ 1 cxp( _ w 2 A 2 ) l x(A)dAj (ig3) 



where 



X (A) = Re ((-a 2 (A)^ Jrfe) - ib 2 (X)\ 2 J[{p 2 )) b* 2 (\)\? J*(p 2 
-Re (^ a2 (A)^ Ji(p 2 ) (-m*(A)A|J(*(p 2 ) - 6|(A)^ J^)) 

= J 1 (p 2 )Jr(p 2 )Re f-a 2 (A)6*(A)^Af + a 2 (X)b* 2 (X)^ 
\ ea e a 

+Re (^J 1 (p 2 )J 1 *(P2)(-^2(A)6;(A)A 2 A; 2 ) + Ji(p 2 )J[*( P2 ) ^a 2 (A)a^(A)^A 2 
Ji(p 2 )Jr(p 2 )Re (a 2 (A)&*(A)^-(-A* 2 + A 2 ) + (^ 2 (A)&*(A)A 2 A* 2 - ia ;(A)a 2 (A)^ 



As e — A 2 = 7 2 , we have A 2 — A 2 2 = e — e*. The following therefore, 
property of the Bessel functions 



yields 



x J[{x)_ x Mx) 
Ji{x) Ji{x) 



(195) 



(196) 



(194) 



X (A) = - a \ J!( P2 )| 2 Re (a 2 (A)6^(A) 7 (l - jj i fajfifi - l) fc* - 7 2 ) (l^(A)| 2 + ^|a 2 (A)| 2 )) . (197) 

I 



There are some reasons to believe that a similar pro- 
cedure can yield an exact formula for the power in the 
incident Gaussian beam, but in this work the power in 
the beam is evaluated using an asymptotic formula for 
Wq ^> 1. As Eqs. (|176H178p define a solution of the free 
Maxwell equation that is accurately approximated by a 
Gaussian beam for wq ^> 1, we shall use the expres- 
sions for the electric and magnetic fields of the Gaussian 



beam at z = (Eqs. (7-8) of (Ref. [1]). The averaged 
p-component of the Poynting vector at the point{a;, y, z} 
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(z — 0) equals 



1 c 



- —Re [E mc {x, y, 0} x H inc *{x, y, 0}] z 



1 r 

T> ^ ( Tpinc Tjinc* rpinc Trine*} 

= 2 47r e v ~ v x ' 

1 c / -\- 

- — Re(i ■ (-i) + 1) exp — 

2 in \ Wn 



47 6XP 



x 2 + y 2 



(198) 



Therefore, the power in the incident beam equals 

x 2 + y 2 



c 
4tt 



oo />oo 



exp 



oo J — oo 
■2tt 



dxdy 



dip / rdr exp 
4tt Jo Jo 



. °° Wndu , , CMJn 

— 2tt / -2— exp(-u) = — ^ (199) where 
47r 7 2 4 

I 



(in Eqs.(7-8) of (Ref. [Tj we replace 5 with wq). 



Thus, Eqs. (|193)197|199[1 yield an expression for heat- 
ing efficiency: 



-21TW 2 . / exp(-wgA 2 )^V(A)dA, (200) 



^A) = |Ji(p 2 )| 2 Re U 2 (A)6|(A) 7 1 



r 



(201) 



Let us find a relation between this expression and Eq. (TTJ 
for wq ;§> 1. The only area that contributes to the integral 
of Eq. (|200[) is that where w X < 1, so in this case we 
may assume that A <C 1 and, therefore, 7 « 1. If we also 
adopt the assumptions of Eqs. (107-110,132) of (Ref. Q, 
we may use the following approximations (see Eq.(116) 
of (Ref. d 



a 2 (A) « 
62(A) 



piJi(p2)piD 

% 



P 2 2 Ji{P2)PiD 



[a + i(3)ie 
(a + 0) 



2 1 



piJi{p2)piD 



IS, 



2 1 



P 2 2 Ji{P2)PiD' 



(202) 



Therefore, 



^(A) « I J 1 (p 2 )n6 2 (A)| 2 Re I ie I 1 - - ) - i I p 2 ^4 " 1 ) fc* " 1)0- + e) 



Vi(p : 



2 



4 D /. ■ * ■ Jo(P2) 1 * . * 1 \ , ■ * • , • * 

1 4 21 nl2 ^ e 1 ?£ ~ Z£ — IP2 7 \ (£ — £ + ££ —!) + «£ — l£ + Z££ 



b 2 | 4 P 2 |C| 2 V "Jifo) 



I 

Let us prove that In fact, as 



4 Im ( P2 ^p\{e* -£ + ££*-!)). (2( ):-!) 



« 2 



rim p 2 — -(£ -£ + ££ - 1) 1 - 



2M 4 V Ji(P2) / (see Eq. © of this work and Eq.(llO) of (Ref. [H), we 

to (± + £±i_iwy ,,04) 

\p 2 2 P2J1CP2)/ 
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obtain that the left-hand side of Eq. (|204[) equals: 

■ 2 ' 1 MP2) 



-Im . 

2 Vpf PiJM 



-Im 



1 



Jo(P2 



(e* -l)a 2 p 2 J!(p 2 ) 



(e*-l)(l + e) 



/!+£ Jofe) 
I +e / 1 . ^(P2) 



= Im 



Im 



2 VpI P2MP2) 

g-1 1 1 1+e J^pg) 

2 (e-l)a 2 p 2 2 ^j^) 

\P2 2 P2JliP2) 



4 2 /l 1 + e J((p2) 
r lm ' 



p?|£>| a a 2 VpI 2 p 2 J x (p 2 ) 



This proves Eq. ([203]) . Thus, Eq. PM)) yields 

^(A) 

If we denote 

m 16 / 1 1 + e J[( P2 ) 

Pi\D\ \P2 2 p 2 Jl{p2) 

then Eqs.(129,143) of (Ref. [1]) yield 

4 



77(A) ~ Im 



therefore, 



^(A) « 17(A) 



(-16)a 2 



-2l(A)A 2 , 



and 



77neu) ~ TTWq / ex P(~ w ^T^M ^ 



Tiwlj cxp(- Wo 2 A 2 )7 ? (A)AdA = |r], (211) 



where 77 is a weighted average of /7(A): 

J °° r?(A) exp(— w 2 A 2 )AdA 
J °° exp(— WpA 2 )AdA 



(206) 



(207) 



(208) 



(209) 



(210) 



(212) 



IV. EFFICIENT HEATING IN THE 
TRANSVERSE GEOMETRY (MORE GENERAL 
CONDITIONS) 

Formulae for heating efficiency in the transverse geom- 
etry were derived in the first part of this work (Ref. [l[, 



Eqs.(156,158)). The formulae were derived in the as- 
sumption that e" ^> |e'|. Let us remove this condition 
now and obtain somewhat more general formulae. Let 
us rewrite the formula for heating efficiency of (Ref. [l[, 
Eqs.(151)): 



\b 2 \R 



4i 



where 



Ai 
naQ ' 



(213) 



(214) 



Q = V^J (V^a)H ( L > (a) - J (^a)H^'(a), (215) 

R = jQ(y/ea)y/eJ (^a)-J {Vea)^/e*J *(y/ea). (216) 
We still have 



R 



-is a 



(217) 



(Ref. [lj], Eqs.(155)), but we need more precise approxi- 
mations of the cylindrical functions of small argument to 
calculate Q: 



Jo(z) 

•to 



# (1) w 



z 2i 
1 

2 7TZ 



2i 

1 + — 



IZ 
7T 



In 



i_f! 
4 ' 

z 

~r 



(218) 



where C — 0.5772 ... is the Euler constant. We leave 
more terms in the expansion of the first derivative of 
the Hankel function as it's magnitude is large for small 
argument. Then, neglecting some terms of higher order 
in a and p — ^fea, we obtain: 



Q 



^ ea 2 \ / a 2i ia^ 
4 / V 2 ira it 



■f(l + |( ta g 1+ C 



a 2i 
2 na 



ia iea ea ( 2% , ,\ 

i+ — = -— 1 + - + 

TT 2n 2 \ tt J 

2i ia^ ea * ea ^ a 
2 ira 7r 2 7r 2 

2i ia a , ia _ , 2i 

L = -- e - 1 + -i(e - 1) = 

7ra 7r 2 7r 7ra 

, (\ L\ 2i , 
= -(e-l)o --»- -— ,(219) 
V 2 7r / 7ra 
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where 



(220) 



Let us find x\ and X2 > that yield the minimum of A 
(and, therefore, the maximum of 77). Differentiating A 
with respect to x\, we obtain: 



Thus, 

U — 4i 

4 

^(-( e -l) a 2(l-ii)-£) - 

4 

~ 2 - i(e - l)a 2 (I - iZ) ~ 
4 

~ 2- (e- l)a 2 (T+if) ~~ 
4 

~ 2 - (a:i + ix 2 ) (L + if) "~ 
4 

2 — Lcei — i_Z/2!2 — i£a;i + 
4 

" (2 - Lxi + §£2) - i (Lar 2 + fail) 
where xi , X2 are real and 

xi + 1x2 = (e — l)a 2 . 



(221) 



(222) 



so 



and 



-!- ^2(2-£a;i)(-L) + ^-a:i ) = 0. 



4L + 2L .ti + —xx = 0, 



21 



Xx 



L 2 



2 

I' 



Differentiating A with respect to £2, we obtain: 



^[{2-Lxxf 



L 2 



7T 

T 



(226) 



(227) 



(228) 



0, (229) 



Then 



16 



(2 - Lxx + f x 2 ) 2 + (Lx 2 + Ixxf 



16 



(2 - Lxx + jx 2 f + (Lx 2 + § xx) 2 



HI) 
(--) 



4tt 



i ((2 - La* + f x 2 ) 2 + (Lx 2 + ^xx) 



(223) 



(224) 



Let us denote the denominator of the latter expression 
by A. Then 



therefore 



x 2 L 2 



so 
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4Lxi 



x\ 



T 2 77 

T 



4 - 4i.xi + 2L.xi = 4 
4L 2 



L 2 



- 2Lx x 

IT 2 

L 2 + 2 



(230) 



X 2 = 



L 2 



(231) 



A = — (2 - L.T1) 2 + tt(2 - Lxi) + ^-x 2 

X2 4 

r2 t 1 71-2 2 

x 2 + irLxx H -a;-, = 

£2 4 



1 

X2 



(2 -in) 2 



i-2;-r- I /. 2 • — ) x 2 . (225) 



ft is easy to see that in the maximum r\ = 1. Reflec- 
tionless absorption of cylindrical electromagnetic waves 
incident on uniform cylinders with small electrical dimen- 
sions was studied in [2j, but the relevant case was not 
mentioned there, perhaps because that article used the 
approximation of Hankel functions for small argument 
that was not applicable in the case of zero index. 

Let us rewrite Eq. (|225|) as follows: 
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A = 




(232) 



For 77 > -J-, we should have or 
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47T 

< T (233) 
I 



,11 2L \ 1 7T 2 . 4tt 

max - xi - — -5- , — JTI' 3 ^ > < -p, (234) 



or simultaneously 



•i'2 




X 2 ( L 2 + 



2L 



L2 + ^, 
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4tt 






T' 
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4tt 
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47T 









(235) 



Finally, after some transformations, we obtain the fol- 
lowing conditions of efficient heating in the transverse 
geometry (these conditions are sufficient, not necessary, 
and describe a domain of parameters where the diameter 
of the cylinder is much smaller than the wavelength): 



x 2 



•t'2 



< 



> 



477 

T' 
1 



2L 



Xi 



< 



L 2 



47T 



X2, 



where x±, xi are real and 

X\ + 1x2 = (e — l)a 2 , 



(236) 



(237) 



a -C 1 is the radius of the cylinder (we consider a 
monochromatic beam with frequency u> and use such 
a system of units that the wave vector in free space 



kn — — — 1); the fields in the incident beam are described 
by (Ref. [1(, Eqs.(144)), e is the complex permittivity of 
the cylinder, 



L 



C = 0.5772. . 



i- h (5 

is the Euler constant 



-c, 



(238) 



V. EFFICIENT LASER HEATING OF 
NANOTUBES 

Let us discuss laser heating in the transverse geome- 
try of carbon nanotubes (e.g. multi-walled nanotubes) 
as the extreme case of ultrathin targets. The resulting 
hot dense plasma may become a lasing medium for gen- 
eration of coherent extreme ultraviolet (EUV) and soft 
X-ray radiation. In our example, lasing may occur as a 
result of recombination in H-like carbon rapidly cooled 
by hydrodynamic expansion (Ref. Q). Laser irradiation 
of multiple nanotubes was previously used for X-ray gen- 
eration (Ref. [![). Let us show, however, that efficient 
laser heating of a single nanotube is also possible. 

Let us assume that the wavelength of the pumping 
radiation is 800 nm (Ti:Sapphire laser) and the pulse 
length is about 30 femtoseconds (3-10~ 14 s), the diameter 
of the nanotube is 15 nm and its density is 1.3 sj cm 2 . 
Nanotubes can be several centimeters long (Refs. @|, 0), 
but, due to the cylindrical symmetry of the problem, the 
following calculations are performed per one centimeter 
of the length of the nanotube. 
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Nanotubes may have very complex electric properties, 
but their atoms will be ionized rapidly in a strong elec- 
tromagnetic field, so we assume that they have electric 
properties of plasma cylinders with complex permittivity 



£ = 1- 



where 



lu(uj + iv ei ) • 

1 

( Aim e e 2 \ 2 



(239) 



(240) 



is the plasma frequency, n e is the electron density (it is 
assumed that carbon atoms are totally ionized), 



Zu p 9N D 

u P 4 ~ — In 

ei 11N D 2Z 

is the electron- ion collision frequency (Ref. [H), 

47T 



N D = Y X D n e 



(241) 



(242) 



is the Debye number, 



(Ref. 0) and resonance absorption (Ref. (sj) should not 
play a major role. 

The energy required for total ionizing and isochoric 
heating to T e = 1 kcV equals 1(1 + Z)T e + Uu per atom 
(although the ion temperature may be lower than the 
electron temperature) , where Uu ~ 1 keV (Ref. 0) is 
the total ionization potential of a carbon atom. Thus, 
the energy required to heat a 1 cm long carbon nanotube 
with a 15 nm diameter to this temperature is about 0.2 
mJ. 

The heating efficiency at T e = 1 keV is about 16%. 
Thus, the required energy in the laser pulse may be es- 
timated as 1.3 mJ, and the laser power is about 45 GW. 
The power of modern Ti:Sapphire pulsed lasers may be 
as high as 100 TW. It should be noted that the heating 
efficiency is higher at lower temperatures. For example, 
the efficiency is about 40% at T e = 200 eV. The effi- 
ciency is about 1% at T e = 10 keV. It is possible to 
achieve higher efficiency at this higher temperature by 
varying the wavelength of the heating radiation and the 
nanotube diameter. It is also possible to use hydrody- 
namic expansion of the plasma. 

VI. CONCLUSIONS 



( T e 



(243) 



is the Debye radius, Z = 6 is the atomic number of car- 
bon, T e is the plasma electron temperature. If T e =l keV, 
the ion sound velocity is of the order of 10 7 cm/s, so the 
radius of the plasma cylinder should not change dramat- 
ically over the laser pulse length. 

The above simple model of plasma electromagnetic 
properties does not take into account nonlocality of the 
complex permittivity. This is partially justified by the 
fact that the only nonvanishing component of the elec- 
tric field is directed along the axis of the cylinder, so 
the normal component of the electric field on the sur- 
face of the cylinder vanishes. Therefore, vacuum heating 



This Part II concludes the study (Ref. 1]) of condi- 
tions of efficient heating of a thin conducting cylinder 
by a broad electromagnetic beam in the longitudinal and 
transverse geometry. 

As a spectacular application, efficient heating of a nan- 
otube by a femtosecond laser pulse is discussed. 
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